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We study the phase diagram of a frustrated spin- 1/2 ferromagnetic chain with anisotropic ex- 
change interactions in an external magnetic field, using the density matrix renormahzation group 
method. We show that an easy-axis anisotropy enhances the tendency towards multimagnon bound 
states, while an easy-plane anisotropy favors chirally ordered phases. In particular, a moderate 
easy-plane anisotropy gives rise to a quantum phase transition at intermediate magnetization. We 
argue that this transition is related to the finite-field phase transition experimentally observed in 
the spin-1/2 compound LiCuV04. 



I. INTRODUCTION 

The interplay of frustration and quantum fluctuations 
in reduced dimensions often leads to unconventional mag- 
netic order, such as chiral or spin- nematic states (see, 
e.g., Refs. HiiiiQIlMIiillillilllilil) • A par- 
ticularly simple model yet realizing a fascinating vari- 
ety of competing phases is the frustrated ferromagnetic 
spin-1/2 chain in the presence of an external magnetic 
field ^i^i^ '^^i^^ described by the Hamiltonian 

= ^ ^Ji{Si ■ Si+i)a + -hiSi ■ Si+2)a - hSi^, 



{Si ■ S2)a = S^S^ + SfSl + AS'iS^, 



(1) 



where Si is a spin-i operator acting at site ^, Ji < 
and J2 > are the nearest and next-nearest neighbor 
exchange constants, h is the external magnetic field, and 
A is the exchange anisotropy. The system may be alter- 
natively viewed as two antiferromagnetic chains coupled 
by a ferromagnetic zigzag-type coupling whose strength 
is measured by the frustration parameter 



13 = J1/J2 



The isotropic (A = 1) version of this model has a rich 
magnetic phase diagram exhibiting states with differ- 
ent types of competing unconventional orders i^s^i^iii 
In particular, the vector chirality, being the quantum 
remnant of the classical helical spin order, competes 
with multipolar orders which characterize the pseudo- 
condensate consisting of multimagnon bound states. A 
similar effect has been previously predicted^ and recently 
confirmed numericall y ^^1^^ for the case of the antiferro- 
magnetic frustrated chain with Ji > 0, J2 > 0. 

The vector chirality (spin current) Ki = (Si x 
Si-^-i) can, even in one dimension, exhibit true long- 
range order (LRO), i.e., the asymptotic value Kq = 
lim|„_„/|^c>o Cnin, n') of the chirality correlator 

C«(ri,n') = «<-) 



can be finite. In the presence of an external magnetic 
field or of a finite anisotropy A ^ 1, the rotational SU (2) 
symmetry is broken down to U{1) x Z2, and the vector 
chiral order corresponds to the spontaneous breaking of 
the discrete Z2 (parity) symmetry. At a finite magne- 
tization, the presence of a nonzero vector chirality au- 
tomatically leads to the emergence of scalar chirality^ 
defined as a mixed product of three spins on a triangu- 
lar plaquette. It has been shown recentljsi^ that in the 
underlying electronic system the presence of a scalar chi- 
rality always induces charge currents, leading to orbital 
antiferromagnetism. 

A common feature of the multipolar phases is that the 
excitations that correspond to a single spin flip (i.e., to 
a change AS^ ±1 of the z-component of the total 
spin) are gapped, and therefore, the in-plane spin cor- 
relator {S^S~,) decays exponentially with the distance 
\n — n'\. This distinguishes such phases from the usual 
spin fluid phases (also called the XYl-type, in the classi- 
fication due to Schula22i) where the {S^S~,) correlations 
decay algebraically. At the same time, the excitations 
with A5^ = ±2 are gapless in the quadrupolar phase, 
those with AS' = ±3 are gapless in the octupolar phase, 
etc. The long-range quadrupolar (nematic) order, char- 
acterized by the finite asymptotic value of the correlator 

C2{n,n')^{S+S++,S~,S-^,) (4) 

at |n — n'\ — > 00, would break the U{1) symmetry, such 
that those correlations can only be gitasi-long-range (i.e., 
exhibiting a power-law decay) in purely one-dimensional 
(ID) systems, yet they may develop into a true LRO 
in real materials where a finite three-dimensional in- 
teraction is always present. The same applies to the 
higher multipolar order parameters such as the octupo- 
lar (triatic) order defined by the correlator of the type 

Cain) = {S^S^+iS^+2Sr+nSr+n+iSr+n+2)^ ^tc. Finally, 
the spin density correlator 

CsDw(n, n') = {S:,S:,) - {S:){S:,) (5) 

has a power-law decay in multipolar phases (as well as in 
(3) the other phases mentioned above) , and depending on the 
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FIG. 1: (Color online) (a) and (b) Magnetization curves for 
P = -0.3 and (a) A = 1.4 and (b) A = 0.5. (c) Phase 
boundaries in the M vs A plane, based on M{h), for /3 = —0.3 
(squares: TV = 96; stars: extrapolation in 1/N). 




FIG. 2: (Color online) (a) and (b) Magnetization curves for 
f3 = -0.6 and (a) A = 0.6 and (b) A = 0.5. (c) Phase 
boundaries in the M vs A plane, based on M{h), for (3 = —0.6 
(squares: A'' = 96; stars: extrapolation in 1/N). 



dominant correlations, a multipolar phase can be further 
characterized as being of the nematic (triatic, etc.) or 
spin-density wave (SDW) type. 

In the isotropic model at small \(3\ the spin gap is pre- 
dicted to be either zero^iiS^ or astronomically smalli^ 
The zero-field phase diagram of the frustrated ferromag- 
netic chain with an anisotropic exchange has been stud- 
ied, both for the case of anisotropic nearest-neighbor in- 
teractions onlji^iiS and for the case in which both cx- 




FIG. 3: (Color online) (a) and (b) Magnetization curves for 
/3 = -3 and (a) A = 1.1 and (b) A = 0.4. The arrow in 
(a) indicates the location where the magnetization starts to 
increase in steps of AS^ — 3. (c) Phase boundaries in the M 
vs A plane, based on M(h), for /3 = —3 (extrapolated in the 
system size TV — » oo). 



change paths exhibit the same anisotropyj ^^i^'''^^i^^ For 
the latter example, relevant to our work, the existence 
of dimer, spin fluid, and (anti)ferromagnetically ordered 
phases has been suggestedi^ Moreover, a chirally ordered 
phase has been predicted to exist at \(3\ < 1 for A < li^ 

The model ([1]) has been suggested to be rele- 
vant for the description of several recently discov- 
ered quasi one-dimensional magnetic materials such 
as LiCuV04 (Ref. M), RbzCuaMogOia (Refjsil), 
Li2ZrCu04 (Ref.m, and anhydrous CuClz (Ref.H^ 

Our goal is to study the interplay between the ex- 
change anisotropy and the magnetic field as reflected 
in the magnetic phase diagram of the model (H]). The 
motivation for our work stems from the experimental 
results^SiMi^^ for LiCuV04 that have revealed the exis- 
tence of a phase transition in a magnetic field from a he- 
lically ordered state at low field values into another phase 
at high magnetic fields where the magnetic order seems 
to be coUinear and directed along the field axisi^^*^ If 
one imagines "switching off" the three-dimensional in- 
teractions, the helical phase might get transformed ei- 
ther into the chirally ordered phase or into a usual spin- 
fluid XY phase (albeit with incommensurate spin corre- 
lations), while the unknown high-field phase could corre- 
spond to the quadrupolar state of the purely ID model. 

However, for the specific parameter values suggested 
to be relevant for this particular material, i.e., [3 « —0.3, 
the ID model ^ with isotropic interactions (A = 1) 
does not support any phase transitions at intermediate 
field values ji^i^ Numerical results^i^ for |/?| > 1 suggest 
that the vector chiral phase shrinks very fast with de- 



3 



creasing \j3\, and thus it is hardly detectable already at 
/3 ~ — 1 . Although one might assume that the vector chi- 
ral phase still persists in an infinitesimally narrow region 
that vanishes asymptotically at /? — » 0, this would not 
sufhce to explain the finite-field transition in LiCuV04 
occuring at a relatively high field strength of about 20% 
of the saturation fieldi^ 

At the same time, electron spin resonance 
experiments^!^ indicate that the exchange inter- 
actions in LiCuV04 have an easy-plane anisotropy 
of about 10%. This puts forward a natural question 
whether including this type of an anisotropy may drive 
the sought-for phase transition. We show that this is 
indeed the case: there is a finite window of A < 1 where 
the spin fiuid phase persists at low fields, while the 
quadrupolar-SDW state occupies the high-field region. 

In the present study, we focus on parameter values rel- 
evant for LiCuV04 (Ref.iOLanhydrous CuClz (Ref.lH), 
and RbaCuzMoaOiz (RefTlHll), namely /3 = -0.3, -0.6 
and —3; respectively. To carry out the numerical analy- 
sis, we employ the density matrix renormalization group 



(DMRG) method 



39.40.41 



and our study is mainly based 



on the calculation of magnetization curves M = M{h) 
and the chiral order parameter kq. In SecHH we present 
the analysis of magnetization curves M = M{h) as a 
function of the exchange anisotropy A. From the magne- 
tization curves, we are able to extract the phase bound- 
aries. The results of the magnetization curves analysis 
are further supported and supplemented by the analysis 
of correlations presented in Sec. IIIII Our main result, 
the magnetic phase diagrams derived from the combined 
analysis of magnetization curves and correlations func- 
tions, is presented in Sec. IIVI We conclude with a sum- 
mary and discussion in Sec. [V] 



II. MAGNETIZATION CURVES 

In this Section, we present magnetization curves of 
the ferromagnetic frustrated chain and discuss their re- 
lation to the phase boundaries of the model ([T|) in the 
magnetization vs anisotropy plane. To that end, we 
compute the ground-state energies Eq{S^) for all values 
S^- = Sf of the z-component of the total spin. By 
subtracting the Zeeman energy —hS^ and carrying out 
the Maxwell construction, we find, for each given field 
h, the quantum number and respectively the magne- 
tization M = 2S^- /N of the ground state, where N is 
the number of sites. Typically, we use about m = 600 
DMRG states, and open boundary conditions. 

Our results for (3 = —0.3, —0.6, and —3 arc shown 
in Figs. [H [21 and [Sj respectively. For both (3 = —0.3 
and P = —0.6, at A = 1 the system is in the quadrupo- 
lar phasei^i^i^ (also called "even-odd" ^ or XY2 phase 
in the classification by Schula^S): the magnetization in- 
creases in steps of AS*^ = 2, due to the presence of two- 
magnon bound states. In this phase, the AS^ = 2 sec- 
tor, corresponding to the simultaneous flip of two spins. 



is gapless, while single-spin excitations with AS^ = 1 are 
gapped.— In terms of correlation functions, at small fields 
the leading instability is in the SDW channel, while at 
high fields the quadrupolar (nematic) correlations of the 
type (HI) dominate 

The results of Figs. [1] and [2] show that an easy-axis 
anisotropy A > 1 simply stabilizes the AS*^ = 2 phase 
[see, e.g., Fig. [Ha)]4^ In contrast to that, an easy- 
plane anisotropy, A < 1, disfavors the formation of two- 
magnon bound states and eventually, we observe the dis- 
appearance of the AS'' = 2 region [see Figs. [Ub) , [D^a) 
and[2l^b)], giving room to the phase with AS^ ~ 1. In 
Sec. mil we will see that this region exhibits chiral or- 
der. The results of the analysis of M{h) for |Ji| < J2 
are summarized in Figs.[TJc) andlDJc): in both cases, be- 
low A < 0.5, the quadrupolar phase has disappeared. It 
is worthwhile to remark that in the case of a weak cou- 
pling {(3 = —0.6 and (3 = —0.3), we observe a reentrant 
behavior in the vicinity of A 0.55: as the magnetiza- 
tion increases, one starts in the AS^' = 1 region, then 
enters into the quadrupolar phase, and reenters into the 
AS*^ = 1 one at M « 0.75. As we shall see below, this 
picture is also supported by the behavior of the chirality 
correlations. 

In the vicinity of the saturation field [M = 1) the 
position of the boundary of the quadrupolar phase is 
in good agreement with the analysis of Ref. |4J. Ac- 
cording to Ref. 0, the field /is,2m at which the two- 
magnon bound state gap closes is given by hg^iml Ji = 
[1 + (A -f- 1)2 - A2(l - /3)2]/[2(l - ^A)] , while the respec- 
tive value for one-magnon states is given by hs^i„-J J2 = 
(A + + + /3)^/8. Comparing those two fields, 

one finds, for example, that for [3 — —0.3 the instabil- 
ity of the fully polarized state at the saturation field 
is by condensation of the two-magnon bound states at 
A > As ~ 0.54, and by one-magnon states below that 
value. The critical point Ag is only slightly dependent on 
the frustration (3, e.g., at /3 = —0.6 one has A., ~ 0.58. 

Let us now turn to the regime of strong coupling, 
(3 — —3. In the isotropic case, the system is in a chi- 
ral phase at small magnetizations, and with increasing 
M one enters a multipolar (actually, octupolar) phased 
This octupolar phase is characterized by AS^ = 3 steps 
in the magnetization curvcj^ which indicates that three- 
magnon bound states are excitations with the lowest en- 
ergy per unit of AS*^. Similar to the < 1 case, an 
easy-axis anisotropy A > 1 stabilizes the AS^ = 3 mul- 
tipolar phase. We illustrate this behavior in Fig. [S^a), 
showing the magnetization curve for A = 1.1. In the 
easy-plane region A < 1, the magnetization curve further 
exhibits a kink-like feature at about AI ^ 0.6, as the ex- 
ample of A = 0.4, plotted in Fig. [3jb), shows. We trace 
this kink back to the incommensurability and the emer- 
gence of multiple Fermi points, following the reasoning 
of Refs. li^lisl . The resulting phase diagram for /? = — 3, 
based on the M{h) curves, is presented in Fig. [3l^c). 

Smirmarizing the results of this section, one can say 
that the main feature, common for all values of the frus- 
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FIG. 4: (Color online) Examples of the chirality correlation 
function {k^{x)k^{0)) at I3 — —0.6 for chains with 128 and 
256 sites, for two fixed points in the (A, M) plane. 



tration parameter /3 considered here, is that an easy- 
plane anisotropy A < 1 gives rise to a mid-field phase 
transition from the AS"^ = 1 "phase" at low fields to a 
multipolar {AS^ > 2) phase at high fields. In Sec. Illli 
we will further focus on characterizing the region with 
gapless triplet excitations {i.e., AS^ = 1) and show that 
it actually contains several different phases. 



III. CORRELATION FUNCTIONS 

In this section, we study the correlation functions, 
complementing the analysis of the magnetization curves 
presented in the previous section. While the multipo- 
lar phases are most easily detected by the appearance of 
the AS"^ > 1 steps in the M{h) curves, the region cor- 
responding to AS^ = 1 can actually contain several dif- 
ferent phases. Indeed, the spin-fluid (XYl) phase of the 
easy-plane spin chain, described by the one-component 
Tomonaga-Luttinger (TLl) liquid, the two-component 
(TL2) spin-fluid phase^^ and the chirally ordered phase 
all have gapless excitations in the AS^ = ±1 chan- 
nel. Thus all those phases will show up as a single 
AS^- — 1 "phase" and cannot be further discerned from 
the M{h) studies. Analyzing the chiral correlation func- 
tion C,i{n,n'), we can identify the chiral phase, and the 
rest of the AS^ = 1 region can be divided into the TLl 
and TL2 phases by the line where a kink occurs in the 
magnetization curve (see Sec. 

Within the multipolar (AS'^ ~ p > 2) phases, an addi- 
tional analysis of correlations is necessary to distinguish 
between the regions with dominant spin density wave 
correlations (SDWp-" phases") and those with dominant 
multipolar (nematic, triatic, etc.) ones. Although the 
transitions between, e.g., SDW2 and nematic is only a 
crossover and not a true phase transition in our purely 
one-dimensional model, such an analysis can be helpful 
in understanding what could be the resulting order in a 
real material, where frustrated chains are coupled by a 
weak three-dimensional interaction. 



A. DMRG methods for the calculation of the 
vector chirality 

Using the finite-size DMRG method^i^ in its matrix- 
product formulation)^ we have studied correlators ([3]), 
(|3]), and ^ in chains of iV = 256 spins. This length has 
been chosen since, on the one hand, it is sufficiently large 
to study the asymptotic long-distance behavior of the 
correlations, and on the other hand, it is small enough to 
ensure that the DMRG calculation converges with a mod- 
erate number m of representative states kept. The typical 
value of m necessary to reach good convergence strongly 
depends on the frustration parameter (3 ~ J1/J2' while 
for (3 — —3, TO = 400 is normally sufficient, at smaller 
coupling {(3 — —0.6) this figure grows to m ~ 600-800, 
and in the regime of weakly coupled chains (/? = —0.3) 
one needs to ~ 800-1200, even for large magnetizations 
M > 0.7 where the convergence is generally faster. The 
correlators ([S])-® have been calculated for a large num- 
ber of ground states in sectors with different . They 
have been averaged over the starting and final positions 
n, n', and contributions with n or n' being closer 
fixed "cutoff" (taken here to be 20 sites) to the chain ends 
have been discarded. Typical chiral correlation functions 
are shown in Fig. [H From such data we have extracted 
the asymptotic value of the correlator which corresponds 
to the square of the chirality Kg . 

A proper finite-size scaling analysis of chirality cor- 
relations is, however, hampered by strong boundary 
effects^ii^'^ that tend to spoil the bulk correlations for 
smaller system sizes. Due to that, it becomes difficult 
to distinguish the chiral LRO from a non-chiral phase in 
those situations where the chiral order parameter kq be- 
comes very small. In such cases, we have complemented 
the finite-size DMRG study with another technique, 
namely, the recently proposed^ matrix-product formula- 
tion of the infinite-size DMRG algorithm (iDMRG) which 
allows to treat systems with finite magnetization (in con- 
trast to the conventional infinite-size DMRG method, 
see, e.g., Ref. [l^. We utilize an algorithm with con- 
served U{1) symmetry to constrain the average magneti- 
zation per unit cell. The convergence rate of iDMRG 
is essentially independent of the size of the unit cell, 
which can be arbitrarily large. The advantage of the 
iDMRG is that the scaling in (m, N) is replaced by the 
scaling in the number of states m alone, which can be 
translated into a scaling with respect to the correlation 
length via ^ to'', where the correlation length ^ is de- 
termined from the next-leading eigenvalue of the transfer 
operator4^ For critical states described by a conformal 
field theory (CFT), is a function of the central charge4^ 
The spectrum of the transfer operator also gives detailed 
information about the exponents and operator content of 
the CFT,^ 

In the standard finite-size DMRG formulation, the de- 
generacy of two chirally ordered ground states will be 
lifted by finite-size corrections. Therefore the purely real 
ground state of a finite system is obtained as a super- 
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FIG. 5: (Color online) Square of the vector chirality Kq: (a) vs magnetization M at (3 — —0.6 and various fixed values of the 
anisotropy A; (b) the same at /3 = —0.3; (c) vs the anisotropy A at two fixed values of A'l and /3 = —0.6 (the arrows here denote 
the respective positions of the AS^ = 1 phase boundary); (d) vs the frustration parameter = — J1/J2 at fixed A = 0.2 and 
M = 0.703125. Open symbols denote the finite-size DMRG results, and solid symbols correspond to iDMRG. 



position of states with zLkq. The IDMRG^^ in contrast, 
allows for a spontaneous breaking of the parity symme- 
try, which also breaks time reversal symmetry and leads 
to a complex valued wavefunction. This gives a transfer 
operator that is not Hermitian, but is instead complex- 
symmetric. The chirality order parameter can then be 
calculated just as kq = 3(S'+5'~^_j). The iDMRG ran- 
domly selects one of the two ground states, with kq ei- 
ther positive or negative. For broken symmetry states 
the iDMRG is quite cfRcient, because the broken symme- 
try state requires fewer basis states than a superposition. 
For example, the representation of a superposition of the 
±Ko states in a form of a matrix product state requires 
precisely double the number of basis states, because the 
reduced density matrices of the two degenerate ground 
states have no overlap in the thermodynamic limit. In a 
finite-size calculation, the mixing of the two states leads 
to somewhat less than a factor 2 in the required basis 
size, nevertheless one still requires generally fewer states 
in iDMRG compared with its finite-size counterpart. 



B. Vector chirality in the weak coupling regime 
(1^1 1 < ^2) 

The results for the weak coupling regime |/3| < 1 are 
shown in Fig. [5] One can see that the finite-size DMRG 
results give the impression that both at /? = —0.6 and 
(3 = —0.3, the vector chiral LRO vanishes in the low- 
field part of the AS^ = 1 region. However, as mentioned 
above, we cannot reliably detect the presence of a very 
small chiral order with the finite-size DMRG method be- 
cause of strong boundary effects. Applying the iDMRG 
technique, one can clearly see that the finite-size DMRG 
tends to underestimate the value of the chiral order pa- 
rameter Kq, cf. Fig. [nila,b)i^ 

As can be seen from Figs.lSJc) and (d). the magnitude 
of the chiral order parameter diminishes quickly when 
decreases, and also when one approaches the boundary 
of the AS^ = 1 region. The convergence in those cases 
becomes very slow. Figure [5] shows the convergence of 
the iDMRG method at a point close to the AS^ = 1 
boundary: a finite chirality is detected when the largest 
intrinsic correlation length ^ of the method exceeds 100 
sites. Taking guidance from the bosonization picture^ii 
it is fair to assume that the chirality can be detected 
only after 1/^ drops below the value corresponding to 
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# states kept (iDMRG) 

FIG. 6: (Color online) Convergence of the iDMRG method at 
the point /3 = -0.6, A = 0.5, M = 0.2. The largest intrinsic 
correlation length ^ of the matrix-product iDMRG method^ 
and square of the chirality kq are shown as functions of the 
number of states kept in the iDMRG calculation. 
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FIG. 7: (Color online) An example of extrapolation in the 
inverse correlation length ^ for the chiral order parameter kq, 
at /3 = —0.3 and M = 0.25, for two anisotropies A = 0.4 and 
A = 0.5 where kq becomes very small, cf. Fig. [5jb). Up to 
m = 1000 steps were kept in those iDMRG calculations. 



the spectral gap in the antisymmetric sectori^ This gap 
becomes very small when one approaches the AS^ = 1 
phase boundary, or the M = linci^ In such cases [see, 
e.g., the low-field region at A > 0.4 for (3 ~ —0.3 in Fig. 
El^b)], one can use an extrapolation in 1/^ to extract the 
chirality kq; Fig. [7] illustrates that this procedure yields 
a finite value of Kq. Continuity arguments suggest that 
tlie entire AS'^ region belongs to the chiral phase, botli 
for P = —0.6 and (3 = —0.3. This is also consistent with 
the theoretical predictions^ of a chiral phase emerging at 
zero field in a wide range of A in the limit |/3| 0, based 
on the analysis of small systems (if the system is in the 
chiral phase already at ft, = 0, it is reasonable to assume 
that the chirality persists at finite field as well). 
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FIG. 8: (Color online) Square of the vector chirality at 
j3 = —3 as a function of the magnetization M at several values 
of the anisotropy A: (a) finite-size DMRG results for a system 
of = 256 spins; (b) comparison of the finite-size DMRG and 
iDMRG results for paths crossing the chiral phase boundary. 



C. Vector chirality in the strong coupling regime 

The behavior of the chiral order parameter in the 
regime of strong coupling (3 = —3, as extracted from the 
finite-size DMRG and iDMRG calculations, is shown in 
Fig. [H] It indicates the existence of a chiral phase that is 
contained inside a relatively narrow stripe 0.5 < A < 1, 
and the rest of the AS*^ = 1 region should belong to a 
non-chiral spin-fluid phase. The presence of a kink in the 
magnetization curves further suggests that this spin-fluid 
phase is in turn divided into the one-component (TLl) 
and two-component (TL2) spin fluid phases, occupying 
the low- and high-field regions, respectively. 

As can be seen from Fig. [8l^a), the transition between 
the chiral phase and the TLl phase is very sharp: the 
chiral order drops from a sizeable value to zero in the 
entire low-field region, when the anisotropy changes from 
A ~ 0.7 to A ~ 0.65. This suggests that the transition 
is of the first order. 

Figure ^h) illustrates that for (3 = -3, the iDMRG 
results agree very well with the finite-size DMRG data 
for a 256-spin chain. This fact, together with the abrupt 
character of the transition from the chiral phase to the 
TLl spin fluid, gives us reasons to conclude that in the 
strong coupling case the observation of non-chiral regions 
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FIG. 9: (Color online) Comparison of the nematic correlator 
C2{x) and the spin density wave correlator Csdw(3;), defined 
by Eqs. Q, ([5}, respectively, as a function of distance x — 
\n — n'l, at /3 = —0.6 and two different magnetizations; (a) 
M = 5/32; (b) M = 25/32. Here, finite-size DMRG results 
for a system of TV = 256 spins are shown. 



is not an artefact of the DMRG convergence, but is due to 
existence of spin- fluid phases, in contrast to the behavior 
in the weak couphng regime |/3| < 1. 
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D. Crossover between the spin density wave and 
nematic at l3 — —0.6 

We have analyzed the crossover between SDW2 and ne- 
matic correlations inside the quadrupoiar AS"^ = 2 phase 
at /? = —0.6. The typical behavior of the SDW and ne- 
matic correlations as defined by Eqs. ([3]), ([S]) is shown in 
Fig. [HI One can sec that both correlators decay as power 
law, but the SDW correlations dominate in the low-field 
region, while the nematic correlations take over at high 
magnetizations, in agreement with the bosonization anal- 
ysis and earlier numerical results j^ii^iii The effect of the 
anisotropy A on this crossover is rather mild: an easy- 
plane anisotropy A < I shifts the crossover boundary 
towards higher M, making the nematic region more nar- 
row, and the crossover boundary seems to be insensitive 
to an easy-axis anisotropy A > 1. 



FIG. 10: (Color online) Phase diagrams of the Hamiltonian 
([T]), represented as slices in the space of the magnetization 
M and anisotropy A, at fixed frustration /? = J1/J2: (a) at 
P = -0.6; (b) P = -0.3; (c) = -3. Solid lines denote 
phase boundaries (extrapolated in 1/N) as extracted from 
M{h). Dashed lines denote the approximate phase bound- 
aries as extracted from the correlations analysis, and dotted 
lines in (a) and (b) denote the position of a crossover be- 
tween the dominant spin-density wave (SDW) and nematic 
correlations. Solid symbols result from the analysis of corre- 
lations functions: Triangles (blue) correspond to the points in 
a chirally ordered phase, diamonds (magenta) to the one- and 
two-component spin-fluid (TLl and TL2, respectively); at 
j3 = —0.3 and (3 = —0.6, circles (red) denote points with dom- 
inant SDW2 correlations, and squares (green) indicate points 
with dominant nematic correlations within the quadrupoiar 
(AS'^ = 2) phase. At /3 = —3, circles (red) denote points in 
the octupolar {AS'' = 3) phase without specifying the domi- 
nant type of correlations (SDW3 or triatic). 
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IV. MAGNETIC PHASE DIAGRAMS 

Summarizing all the information extracted from the 
analysis of magnetization curves and correlations, one 
can establish the phase diagrams of the anisotropic frus- 
trated ferromagnetic spin-i chains in the presence of a 
magnetic field. Such phase diagrams in the {M, A) plane 
at different values of the frustration (3 = Ji / J2 are pre- 
sented in Fig. [ini 

We reiterate here that we ascribe the entire AS*^ = 1 
region to the chiral phase for fi = —0.3 and /3 = —0.6, 
based on the very smooth character of how the order pa- 
rameter vanishes approaching the AS*^ = 1 boundary, on 
theoretical estimates for zero-field case,— and by invok- 
ing continuity arguments. In principle, from our data, 
we cannot exclude the existence of a small non-chiral re- 
gion in the low-field part of the phase diagram near the 
AS*^ = 1 boundary, but we think that this scenario is 
rather unlikely. Thus, in the regime of weakly coupled 
chains, the phase diagram contains just the chiral and 
quadrupolar phases, the transition between them being 
likely a smooth (second-order) one. 

In the strong coupling regime, /3 = —3, our results sug- 
gest a rich phase diagram, displaying four phases: the 
octupolar, the chiral, and two types of spin-fluid phases 
(which can be characterized as one- and two-component 
Tomonaga-Luttinger liquids). The transition between 
the chiral and the spin-fluid-TLl phase is very sharp and 
is likely first order. 

V. SUMMARY 

Motivated by recent experimental results 
for several quasi one-dimcnsional magnetic 
materials,^2i^i22iSi2ii^i2& we studied the model ^ 
of an anisotropic frustrated ferromagnetic spin-i chain 
in an external magnetic field, at finite values of the 
magnetization. 

We showed that an easy-axis anisotropy A > 1 
stabilizes multipolar phases j^i^ii^ in which the total z- 
projcction of the spin S'' increases by steps of AS'^ > 1. 
In the presence of even a small easy-axis anisotropy, such 
phases occupy the entire range of finite magnetizations 
up to full saturation. Further, we found that an easy- 
plane anisotropy A < 1 may favor several types of phases: 
chirally ordered and non-chiral one- and two-component 
spin fluids. We showed that the presence of a moderate 
easy-plane anisotropy leads to the possibility of a field- 
induced quantum phase transition at a substantially large 
value of the magnetization Af , even in the purely one- 
dimensional model IJ]), which might provide an explana- 
tion for the field-induced transitio n'^^'^^i'^^ from a heli- 
cally ordered to a coUinear state observed in LiCuV04. 

Assuming that LiCuV04 is a system of weakly cou- 
pled onc-dimcnsional chains and further assuming that 
the presence of an exchange anisotropy drives the exper- 
imentally observed mid-field phase transition in this ma- 



terial, our results imply that the low-field region would 
be in a helical cone-type phase (see Ref. [H for a recent 
study of helical order in a 3D magnet in high magnetic 
fields) while the SDW instability in the high-field region 
would turn the high field region into a collinear, mag- 
netically ordered state with long-range incommensurate 
{S^{x)S^{0)) correlations. The former conclusion (a heli- 
cal phase in the low-field region) is in agreement with the 
available experimental datapSii^ii^ while the latter con- 
jecture of coUinear incommensurate order could be tested 
by neutron scattering experiments and is consistent with 
the nuclear magnetic resonance measurements^^*^ sug- 
gesting that the magnetic order becomes coUinear in ap- 
plied fields above « 7.5 T. It should be mentioned that a 
similar incommensurate coUinear structure has been re- 
cently observed^ in the quasi-lD material BaCo2V208 
with easy-axis anisotropy. 

For f3 = —0.3, our data suggest that the low-field, chi- 
ral region opens up at a finite anisotropy; within the 
numerical accuracy of our calculations, we were able to 
resolve the emergence of this region for A < 0.8 [see 

Fig.iinrb)]. 

This has to be contrasted against the experimental es- 
timate of the easy-axis anisotropy of about 10%r^^ and 
against the fact that in the magnetization measurements 
in LiCuV04, the mid-field transition is observed at 7.5 T, 
corresponding to about 20% of the saturation fieldi^i^i 
We stress that our results do not serve to unambigu- 
ously prove the exchange anisotropy to be the relevant 
mechanism behind the mid-field transition in LiCuV04; 
nevertheless, our results clearly indicate that, using the 
values for (3 and A suggested for LiCuV04, this material 
is very close to the quantum critical point at which, as 
a function of decreasing A, a mid-field phase transition 
develops. This transition point shifts to larger field as 
the anisotropy increases (A decreases). 

The vicinity to many competing phases then makes 
this material so interesting but also renders it difficult to 
quantitatively predict its phase diagram. Additional ex- 
perimental data are highly desirable to clarify the nature 
of this phase transition, while, in conclusion, our work 
shows that the emergent physics in this model, driven 
by the magnetic field, quantum fluctuations and broken 
exchange symmetry, is very rich. 
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